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METHOD FOR APPROXIMATING THE VACUUM MOTIONS OF SPINNING 
SYMMETRICAL BODIES WITH NONCONSTANT SPIN RATES 

By C. William Martz 


SUMMARY Jtf ^ $/ 

A method for approximating the vacuum motions 
of spinning rigid symmetrical bodies with varying 
spin rates and inertias has been completed. The 
analysis includes the effects of time varying thrust 
misalinements , mass unbalance , and jet damping. 
Results are given in the form of equations for space- 
referenced Euler angles , flight-path angles , body- 
referenced attitude rates , and earth-referenced 
vehicle-trajectory coordinates. The method consists 
of dividing the problem into intervals during which 
the time-dependent variables are assumed constant 
at their mean interval value. In order to check this 
procedure , solutions for various interval sizes are 
compared with solutions obtained with numerical 
methods. Although the method is somewhat lengthy 
for accurate hand computation in most cases , it is 
readily programed for machine solutions. Probably 
more important , the general solutions give insight 
into the separate effects of the variables and , in many 
cases , can be quickly used to determine the approxi- 
mate ranges of the variables required for the desired 
solution to a given problem. In this respect , equa- 
tions for determining maximum wobble have been 
derived for certain input conditions. 

The method has been shown to compare closely 
with the numerical solutions of two sample problems. 
The sample problems also illustrated the relatively 
large effect of pitch and yaw jet damping on body 
motions. 

INTRODUCTION 

Vacuum motions of rotating bodies are becoming 
more important with the fairly recent ability to 
place objects in motion beyond the atmosphere. 
Machine computer programs for calculating these 
type motions have been completed and used 


successfully lor some time. However, not every- 
one has a computer machine available for this 
work. Also, those with machines are using the 
trial-and-error process in most instances when 
locating the proper range of variables with the 
result that much machine time could be saved if 
some insight were available as to the individual 
effects of the different variables on the motions. 
This insight is best provided by analytical solu- 
tions to the equations of motions. There have 
been many papers published concerning this prob- 
lem. (See, for example, refs. 1, 2, and 3.) How- 
ever, one thing common to these papers has been 
the constant spin rate requirement. Other re- 
quirements sometimes include constant mass and 
inertia parameters or constant moment inputs. 
Solutions are sometimes limited to angular rates 
referred to a body-axis system requiring trans- 
formation and numerical integration to obtain 
space-referenced attitude angles. 

The present paper presents an approximation 
method for determining the vacuum motion of 
spinning symmetrical rigid bodies with changing 
spin rates and inertias including the effects of time 
varying thrust misalinement, mass unbalance, and 
jet damping. Results are presented in the form 
of equations for space-referenced Euler angles 
and flight-path angles, and earth-referenced 
vehicle-trajectory coordinates. An expression 
for body-referenced attitude rates is included for 
convenience. The method consists of dividing 
the problem into intervals during which the time- 
dependent variables are assumed constant at their 
mean interval value. In order to check this 
procedure, solutions for various interval sizes are 
compared with solutions obtained with numerical 
methods. The method was developed under the 
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limitations that body pitch and yaw attitudes are 
restricted to “small angle ’ ■ 9 oscillations and that 
body moments of inertia about the pitch and yaw 
axes are equal. 

SYMBOLS 


a 

arbitrary fitting constant 

A, B,C 

complex input coefficients defined in 

A n) Bn 

equation (13) 

fitting constants for moment inputs 

a 4 ,a 5 

constants defined by equations (18) 

B 4} B 5 , 

constants defined by equations (18) 

b 6 ,b 7 

Cu c 2 , 

constants defined by equations (18) 

a, c< 

F 

magnitude of total asymmetrical 

X m 

force on vehicle 

input coefficients defined in equation 

F 

(13) 

mean value of FI {mV) within an 

g 

interval 

constant of gravitational acceleration 

Ixi Iyj Iz 

moments of inertia about X-, F-, 

I XY , Ixz 

and Z-axes, respectively 
products of inertia due to unbalance 


I=Iy=Iz 


x e , Vc, z e x-> y-> anc l z-axis system rotated 
about y-axis to make z e - axis ver- 
tical 

a bod} 7 angle of attack referred to a 

rolling body-axis system 
body angle of attack referred to a 
nonrolling body-axis system 

0 body angle of sideslip referred to a 

rolling body-axis system 

p s body angle of sideslip referred to a 

nonrolling body-axis system 

i 3=j-\-i(p—a>) 

y 0 angle between a?-, y-, and 2 -axis 

system and x e , y e , and z e system 
in £ 2 -plane 

7=7 a +n* 

y e flight-path angle in pitch plane 

7 ^ flight-path angle in yaw plane 

X vehicle total yaw angle, 0+?V, radians 

rj Y) rj z angle between body principal X-axis 

and X (body reference) axis 
measured about F- and Z-axes, 
respectively (see figs. 3 and 4) 

(j> 0 = j * pdt at t 0 


i= V — 1 

j jet damping coefficient, K/I 

J mean value of j over the interval 

K pitch and yaw jet damping factor, 

I -ml 2 


K ' 
l 


m 

M x , M Y , Mz 


P, Q, r 

V 

R„ R 2, R 3 

T 

T 

t 

V 

X,Y,Z 
z, y> z 


roll jet damping factor 
distance from body center of grav- 
ity to motor nozzle exit measured 
along X-axis 
mass of body 

asymmetrical moments about X-, 
F-, and Z-axes, respectively 
angular velocity about X-, 1 -, and 
Z-axes, respectively 
mean value of p within an interval 
vectors defined in equations (9) 
thrust 

mean value of T/(mV) within the 
interval 

time from beginning of interval 
velocity of bod} 7 along flight path 
orthogonal body-axis system (origin 
at body center of gravity) 
orthogonal space-axis system (origin 
arbitrary) 




0' 


Subscripts: 

o 

/ 

max 

n 


yaw, pitch, and roll orientation 
angles of body X-, F-, and Z-axes 
with respect to x, y, and z space- 
axis system (Euler angles) 
angle between the total asymmetrical 
force vector (always in the FZ- 
plane) and the — Z direction (see 

% 1(b)) 

mean value of pl x l I within the 
interval 

value of quantity at beginning of 
interval 

value of quantity at end of interval 
maximum value of quantity 
integer 1, 2, 3 


A dot over a symbol indicates the first derivative 
with respect to time; a double dot indicates the 
second derivative with respect to time. 

ANALYSIS 

The modified Eulerian dynamic equations gov- 
erning the rotational motions of a body about its 
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principal axes are (see refs. 4 and 5) : 

M x =IxV — Tq_{Iy — Iz)-\~K'p ( 1 ) 

M y =Iy<L — rp(I z — I x )-\-Kq (2) 

M z = I z f — pq(Ix — Iv)+Kr (3) 

Figure 1 illustrates the axis system used. 


If the body is assumed to have rotational mass 
symmetry, I z will be equal to I Y and the rolling 
motion will not be affected by the pitching and 
yawing motions. This allows equations (2) and 
(3) to be solved independently of (1) for pre- 
selected p histories. 

By multiplying equation (3) by i and adding 
the result to equation (2) with the rotational 
symmetry assumption, the equation becomes 



(a) Pitch and yaw orientation of axes systems. 0=0. 



(b) Roll orientation of axes systems. 0=0=0. 
Figure 1. — Axes systems employed in analysis. 


M Y +iM z =I(q+ if) + ip ( I— I x ) (q + ir) 

+K(q+ir ) (4) 


This equation can be referred from a rolling 
body-axis system to a space-axis system with the 
transformation equations (ref. 6) 


6 = q cos <t>— r sin </> " 

H^b) (2 siD <t,+rcos - 

<f>=p-\-\// sin 6 J 


( 5 ) 


Now, for small values of 0 when cos 0« 1 and 
\j/ sin 6<^p (zero reference for 0 can be changed 
when necessary), equations (5) result in 


\=d+i\//= (q+irfe** (6) 

where 

Combining equations (4) and (6) yields 
where 

. K I-rhl 2 

I~ I 

Equation (7) then governs the pitching and 
yawing motions of rotation ally symmetric bodies 
referred to a space-axis system. The general form 


of solution for this equation is 


\=-6+it 




dt-]-\ 0 


( 8 ) 
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The problem now is to find time functions for 

the variables (j , p, w j^ c ] 1 not on ly 

permit equation (8) to be evaluated but which also 
adequately approximate the time histories of 
these variables as the}" would exist in any practi- 
cal problem. 

SPECIFIC SOLUTIONS 


Solution with variables constant. — The solution 
I x M y Mz 

of equation (7) when J,p,-j> ~j~> and are con- 
stants is 


X — 6 -T i\p 

=R 1 +R 2 e i dT)‘+R 3e <*« (9) 

where 

X 0 (M Y -\-%M z )e l<t>0 / 1 

I 3 \ I 

(9a) 



\ 0 e-» (Mr+iM z )e'*°- 1 ‘ / 1 \ (C)h) 

2 ipl x T p 2 (I—I x ) — Upj Ux_j) 

1 J \I ipj 


_ — (M y 4- iM z ) e < *° 
3 p\I-I x )-Upj 


(9c) 



This X solution can be thought of as the sum of 
three vectors: a nonrotating trim vector R,, a 

vector R> rotating at the rate and a vector 


R 3 rotating at the rate p. This type of motion is 
referred to as “tricyclic” in reference 1 and illus- 
trated in figure 2. The low-frequency vector is 
called the precession vector, and the high-fre- 
quency vector is called the nutation vector. Note 
that jet damping attenuates only the R_> vector. 
Equation (9) may be more familiar with j= 0 and 
with the real parts separated from the imaginary 


parts as follows: 


e=e 0 +- e -f- sin p 1x1 ' *• 


v- 


I ' lx 

p 7 


/ Ixt X 

\ , Mr cos <t> 0 —M z sin 4>„ 1" 

1 COS p — J — I 

> 1 pHI-lx) L 


-cos pt 




, It Ixt ,YI M z cos 4> 0 -{-Mr sin <t>„ T . . | I ( 

+rA coap T~ 1 )i W=i3 — 

f— siup(+^ 


sm p 


tl, 


)] (10a) 



\ i i'o ■ Ixt . 

M y cos <t) 0 —M z sin 4> 0 1" 

(cos p j 1 ^ 

J 

K«| 

J 

p 2 (I-Ix) [ 





r) -^Y] 

M z cos <t> 0 -\-Mr sin 4>„ 


\ ln P I )j 

P 2 (I~Ix) 


Exact solution with nonconstant spin rates.— Of the many attempts to satisfy equation (8) by 
substitution of various time functions for the variables, the one which permitted an exact solution 
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with nonconstant spin rates used the substitutions 

Vo 


1 1 +at 

Ho 

lx 

-j-=Constant 


( 11 ) 


The assumed straiglit line dependence of 1/p is 
quite practical especially if the problem is divided 
into time intervals. Although the constant-inertia- 
ratio requirement can be circumvented (by the 
method of the next section), the zero jet damping 
limitation is considered serious except, of course, 
for the case of nonthrusting vehicles. 

Because a more general solution (presented in 
the next section) was found for equation (8) and 
in order to reduce confusion, the exact solution 
referred to in the present section is presented in 
appendix A along with all further discussion of 
this solution. 

Mean value solution. — In the application of this 
solution, the problem is first divided into time 
intervals. The number and duration of these 
intervals depends upon the accuracy desired and 
will be discussed in the section entitled “Results 
and Discussion.” Within each of the time inter- 
vals, the variables p, pl x !l, and j are approxi- 
mated by their mean value over that interval. 

For example, consider the damping term in the 

exponentials of equation (8), namely, /; jdt. 

This integral is approximated by jt where ] is the 
mean value of j over the interval. By definition, this 
is an exact approximation when the integration ex- 
tends over the complete interval ^ i.e., f jdt =jt\ 

For times less than one complete interval, however, 
the result is approximate. The accuracy of this 
approximation can be increased to any desired 
level by using shorter time intervals. Thus, with 
the substitutions 


p=p 

v'-f=z 


3=3 


( 12 ) 


a straightforward integration of the exponentials 


of equation (8) can be accomplished. 

Concerning the moment inputs of equations 
(7) and (8), M Y and M z should be approximated 
by functions which can adequately describe the 
variations of known time-dependent moment 
asymmetries such as thrust asymmetries, tip-off 
asymmetries, and dynamic unbalance effects. 
Remembering inertia must also be allowed to 
vary with time, the following input forms are 
assumed for each interval: 

My+iMz = ^ [(F n +iOn) ( 1 +AJ +B n t*)\ 

=A+Bt+Ct 2 (13) 

When thrust or tip-off asymmetries are con- 
sidered, M y and M z are the actual pitching- and 
yawing-moment asymmetries applied to the 
vehicle. 

When dynamic unbalance effects are considered, 
the moments M y and M z are related to the perti- 
nent variables as follows (ref. 4): 

M Y + iM z = I xz (f 2 p 2 ) 4" I x r ( 7 * T <z) 

+i[Ixz(p— qr ) +Ixy(p 2 — q 2 ) ] (14a) 

which, for the present purposes, reduces to 


My+iM z — p 2 ( IxZ-\-Uxy) ~p 2 (I — lx ) ( VY~\~iVz ) 

(14b) 


since the products of inertia are related to angular 
deviations of the principal axes as follows: 


and 


tan 2rj Y 


tan 2 rj z 


27 xz 
lx Iz 


2/y 


Now, if a combination of asymmetries and un- 
balance exist during the same interval, it may be 
easier to fit each asymmetry or unbalance to a 
separate complex input term. For example, the 


unbalance input p 2 y y— 1 j(rj Y +ir} z ) may 


fitted to the term {Fx+iG^il+Axt+Bffi. If the 
input moments have large or rapid changes in 
direction during an interval, however, it is more 
satisfactory to combine the real components of 
the various inputs separately from their imaginary 
components. Then, the total complex input is 
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fitted to a combination of two or more input terms 
as 

total = 0^1 + Oi) (1 -\-A\t-\- Bit 2 ) 

-b (O+iGQ (1 ~ b 

Now, when equations (8), (12), and (13) are 
combined, equation (8) becomes 


\z= e f *» ‘ CA+Bt+Ct*)e* t+i{ *~ Zu dt 

+X»e-**oJ dt+X, (15) 

By integration, 


p i(pt + <t> 0 ) 

*-r 




(16a) 


and 




ipft 

+ 


where 


P=j+i(p— «) 


(16b) 


Equations (16) predict the approximate rotational 
motions (space referenced) of spinning sym- 
metrical bodies with changing spin rates and 
inertias including the effects of time varying 
thrust misalinements, mass unbalance, and jet 
damping. As in equation (9), the X solution is 
tricyclic with a nonrotating trim vector, a vector 
rotating at the rate oj, and a vector rotating at 
the mean spin rate ~p. By separating the real 
and imaginary parts of equation (16b), thesolution 


mav be more readily evaluated in the following form: 


cos W <- 1 ) (C, w +Cj)( g y> sin a>n + , v , cos pt-U-Ci) sin pt+C> 

co 2 +j 2 ^ 2_ b j 2 

+si [A COS (pt+^+A, sin (pt+<f>o)] f A n t+Bj*- =r^^— l-[A sin (pt+<t> 0 ) 

72 = 1 h L J 

_ , „f 2 B n (v-Z)t . 2 B„n\ 

-A, co, <?<+*.)] 


(17 a) 


. (-C ] ]+C 2 w)(e- i ‘ sin cop (C^+C3)(e » cos *»t-i) + (a gin pt+ih-C*) cos pt+G* 

&+f " 2 +j 2 

+S -f [^4 COS (pt+<t> 0 )+A t sin ~\a[A sin (jpt+<t> 0 ) 

«=i l. U +(p— ») V J 

-A, cos ( pt +4>„) ] ^A n t + B„t 2 - ( 17b ) 

where 

A Fn(p u)-\rG n j (18a) 

Plj 2 +(P — ^) 2 1 

a = f 3+ G *(P— 3. (18b) 

Mi 2 +(p— “) 2 1 

c x =e„+i: (-£ 4 b 6 +5 6 b 7 ) 


(18c) 
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C-2— 


(18d) 


r ff ,4yf-Bv f p 2B„ 2($-u)B n I 

_£ a r /? , A, 2B n ~j -p 
PL 5 P P?+P(P-^)UJ 

c.=*. + £{^lB. + f- W- 1 

»=i L p L p w 2 +y(p— «) 2 J 

■ B, r B _ 2 B„ _ 2 BJ^-w) 1 \ 

p L 4 p 2 pj 2 +p(p-zyJJ 

AJ , 2B,\l i -(,p-«y 


(18e) 


(180 


B<= 1- 


i 2 + (p_-)2 [j2 _ ( p_ S ; )2] 2 +[2i( p_- )]2 

(18g) 


B _ Ar,{p—Z>) 4 . jB„(_l>—u) _ 

5 J 2 +(7>-«) 2 t?-(p-5)T+[2j(p-5)] 2 

(ISh) 

B^J>(A, cos <*>„— A, sin 4>„) (18i) 

B 7 =p(A 5 sin fo+Ai cos <#>„) (18j) 

Because body motions referred to a body-axis 
system are sometimes desired, equations (16a) 
and (6) are combined to obtain lor reference 

q+ir^{q u +ir 0 )-l )J 

or 

g+ir= ( 1 9l>) 

FLIGHT-PATH ATTITUDE 

Up to this point, all effort lias been toward 
determining the attitude of the body expressed 
by the Euler angles 6 and \p. Of greater impor- 
tance to many investigations is the knowledge 
of how the external forces and moments on the 
body affect its velocity vector and space location. 
The author’s interest in an analytical solution 
to this phase of the general problem was stimulated 
by the analytical results of reference 3. The 
method of reference 3 will now be used to extend 
the attitude solution of the present paper to 
expressions defining the velocity vector. 

The force equation normal to the flight path 


in the pitch plane (fig. 3) is 
mVy =T sin a,— mg cos (y+y'„) 

v 0 

T I 1 cos ol s cos (20a) 

and in the yaw plane (fig. 4) is 


mVy^ =—T sin p s +F cos (3 S sin 

(20b) 


By restricting the results to angular changes 
in velocity vector due to external disturbances 
other than gravity, the weight term of equation 
(20a) can be dropped. Multiplying equation 
(20b) by i, adding the result to equation (20a), 
and considering only small angles 3 uelds the 
following: 

(y,+i'Y+) + (7e+iy+) ~y^~y ( 6 + i 'P) 


where 


_Z 1 

mV 



pdt + 0 


■) 


( 21 ) 


0=«»+ 7« 




Figure 4. — Yaw-plane forces and angles. 0 =^= 0 . 
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and 


^ = — &+7* 


The form of solution for equation (21) is 

r 


J 't T 

0 mV 





I 9 ***'**) (Jo' 



dt-\-\ 0 


Again, the mean value substitutions within the intervals are used. Let 


-V* 

mV 

mV^ 


( 22 ) 

(23) 

(24) 


Substituting equations (16), (23), and (24) into equation ( 22 ) and integrating gives 

TV*. r - - r, B.2C, 2 C ,/„ 2C,2iC\/i 1 \"1 

^ip0(T+ip)\ L /S j8* (l+ip)* \ P V J\P 1 +wJ 

"■ . B,2C\ 

l( I + lO 


\ P p T+ipJ J 


10)— J 


[e~ 7 ‘-l _ (,'“'-"''-e- 7 ')] 

7— J+7W J 




+lr ^ -^-"1 { x --S[ / 1 _ f + (ri)( B_ T)]} (25 


5) 


Equation (25) predicts the direction ol the ve- 
locity vector for vehicles having the angular 
motions described by equations (16). Note the 
similarity of these two solutions (eqs. (16b) and 
(25)). Both consist of a fixed vector, a vector 
rotating at the mean roll rate, and a vector ro- 
tating at the mean value of 

Again, it should be remembered that these 
solutions are for small values of 6 and that flight- 
path curvature due to gravity is disregarded. 
An approximate change in flight-path angle due to 
gravity is 

(ft COS 7 o -rt 


and because of the relationship 

f5-\-ia= (0 s ~\~ioi s )e l4> 
a= — / 3 . s - sin (/?/+<£„)+«* cos (pt+<t>o) 
P=& cos (pt+<l>o)+ot a sin (pt+<t>o) 


(27) 


A(7«+nJ=“ 


V 0 


The results of equations (16) and (25) can be 
combined to yield time histories of angle of attack 
and angle of sideslip as follows: 


ot s =d—y e I 
&=— 1^+7* J 


(26) 


Now, going on to the space-position solution, it 
can be shown from figure 1 that 

z e =—V sin ( y e +y'o ) 

x e =V cos (y d +y'o) cos 7 ^ V (28) 

if e = V cos ( t 0 + 7 o ) sin 7 ^ , 

The foice equation along the flight path for small 
disturbances is 

T cos a cos $—m(f sin (y d +y o) = mV (29) 

By considering small angles lor a and /3 and re- 
moving gravity effects equation (29) becomes 

T = V 
mV V 


(30) 
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Integrating this equation and combining with 
equation (23) results in 


V=V* T < 


(31) 


Now, this velocity expression can be substituted 
into equations (28) with the assumption that y e 
and y+ are small angles, and the equations can be 
expanded. Since equations (28) and (31) were 
obtained by neglecting the force of gravity, the 
term gt is added to the equation for z e to get the 
earth-referenced velocity equations with gravity 
effects included. 


z e = — V 0 e Tl (y e cos 7*+sin y' 0 ) + (jt' 
x e =V 0 e Tt (cos y' 0 — y e sin 7;) ► 

Ve=V 0 e Tt y yl/ (cos y ' 0 — y d sin y' 0 ) 


(32) 


\ 0 and the initial attitude rate \ 0 can readily be 
demonstrated by making A=B=C= 0 and ]= 0 
with the result that 


\=X»+4 — 1) 

ICO 

This solution is shown in figure 5(a) to be the sum 
of a fixed vector and a rotating vector. Until 
\ 0 and \ 0 are specified, however, it is not possible 
to say whether a higher or lower spin rate or inertia 
ratio will increase or decrease X max , the common 
performance standard of spin stabilization. 

Now, if the model is thrusting but with no 
disturbance moments, jet damping normally 
attenuates the motion and the plot of 6 against \p 
turns into a logarithmic spiral as shown in figure 
5(b). 


Equations (32) are integrable and yield the 
space location equations with gravity effects in- 
cluded. The result is given here as functions of 
y e and 7^ which are available from equation (25). 


z e =z e ,- V “ S " 1 y '° (e Tt -\) 


— V„ cos 7 oj y e e T ‘dt+\ gt 2 


x e =x e , 0 + v ° c ° sy '° (e r '-l) 


— F 0 siii7' (7 ,e T ‘dt 

=ye. 0 +V„ cos 7 ' J 7 +e T, dt 

— V„ sin 7 oj y,y^e T ‘dt 


(33) 


RESULTS AND DISCUSSION 

The results of this paper are primarily the 
attitude solutions expressed by equations (16), 
the flight-path direction expressed by equation 
(25), and the space-location solution of equations 
(33). All these solutions are complicated by the 
large number of variables which affect the end 
results. In order to show some of the more im- 
portant interrelationships of these variables, 
others must be held constant. 


ATTITUDE SOLUTION 

With the attitude solution of equations (16), 
the separate motion effects of the initial attitude 


9 




(а) j= 0 (no jet damping). 

(б) i>0 (with jet damping). 

Figure 5. — Sample attitude solution with no moment 
inputs (F n -\-iG n = 0). 
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Moment asymmetries. — For the effects oi mo- 
ment asymmetries the solution for X may be 
simplified by assuming no residual motion (\ 0 
= A w =0) and constant inputs (/i-^G r =0). Then 
for <t> 0 — 0, the solution becomes 


x- r^-n ^ 1 /"V'H 

ilp[j+i(p - «)] L w— 3 J 


Expressions for maximum wobble can be obtained 
from this equation under the following conditions : 
When <P, or, more accurately, when 


!=[ i+«‘ 

2 _L_ ^2 


X a) 2 +j : 


-+ --?* - 1 

cu(/J-o>) ^>2 


then 


My 2 +M/ 


^ max 


(u 2 +M 2 +(v-*) 2 } 


^CVc-iy+(|y 

“■+=7^-1 

CO oo(p — co) I 


+X“Si' + / 

V 


ai 1, for normally small j values (i.e., j/p<^.l and 

5/«« D- 


. 2 JM7+M? 

maz ~ Iu(p-<b) 


When co >/; . or, more precisely, when 

Vs 2 +? 

then 

m7+m? 
i 

(«*+?) [?+(p-^ 2 ] 

| V^M-j 2 | e -u/p)tan~Hl/~) J 

and, for normally small j values, 

2 VmT+m? 

" iax ~ /?(«-*) 

Note that both A max expressions show that max- 
imum wobble due to asymmetrical moments is 
proportional to the size of the moment and, for 




normally small J values, is effectively inversely 
proportional to the product of the mean roll in- 
ertia and the square of the mean spin rate over 
the interval. It follows, therefore, that spin-rate 
magnitudes increasing with time (as well as larger 
spin rates) will reduce wobble because of the re- 
sulting larger value of mean spin rate over the 
interval. 

Lastly, for «co, an expression for maximum 
wobble can be obtained under the reasonable 
conditions that (j/p)^ 1 and that ( p — c o)t is a 
“ small angle/ ’ This expression is 


*=iW (Mr ’ +Mz2) 


(l-g-' f, ) 2 +[e- /t (u-J)tY 

j 2 + (p—u ) 2 


This solution indicat es_ the divergent nature of 
A max for p^u J at small jt values. The divergence 
is more apparent for j—0 with the result 


\ max =±jM?+M? 

Although the quantity (p—co) is never exactly 
equal to zero in a practical problem, the theo- 
retical possibility of I=Ix also results in the sim- 
ple divergence equation just given (for normally 
small values of J). These divergence solutions 
reveal that wobble buildup is proportional to the 
input moment disturbances and to the time re- 
quired for passing through the resonant condition 
and inversely proportional to pitch or yaw 
inertia and mean spin rate. 

Unbalance. — It is difficult to show clearly the 
effects of unbalance on the motion of a body hav- 
ing a nonconstant spin rate because the input 
moments (eq. (14b)) are variable. Thus, the 
quantities B and C or A n and B n must take on 
values other than zero which precludes a simpli- 
fied version of the general solution. However, by 

considering fixed p and ^ values temporarily, 

equations (9) and (14b) can be combined (still re- 
taining the conditions A o =A o =0, B = C=0, and 
<t> o =0) with the result 


rjy-{-ir\z 

i iL ! 

rV-ix) 


e ivi — \-\-ip 
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T x 


When -j <1, X maI = 


Vr 2 + nz 2 


V t 1 +(^=fe) ] [(r ) +(?) ] 

and, for normally small values of j, 

21 


(l+etf) 


k max — j V V y 2 I 




When (*>1 


and, for normally small values of j, 


/ Vy 2 +Vz 2 

/ ! + 

r a i 

1 

* 


When ^=1, 


Klax~2^^Y 2 +Vz 2 


These results are for constant values of p and -j 

with no residual motions. Note that wobble due 
to unbalance is independent of spin rate except 
through the jet damping terms which normally 
have only a very small effect. Thus, wobble due 
to unbalance cannot be controlled by spinning as 
can, for example, the wobble resulting from initial 
attitude rates or thrust misalinement. 


APPLICATION 

The analytical expressions for X and y given in 
the present paper have been programed for use 
with an IBM 7090 electronic data processing 
machine. Sample problems were composed and 
this program used to generate their solutions. 
The numerical solutions to these problems were 
also obtained with the numerical integration 
method reported in reference 3 for the purpose of 
comparison with approximate solutions. 

lt_should_be mentioned that the values for p 
oo, T, and F used in the approximate solutions 
were based on assumed exponential histories of 

the variables within the interval I i.e., for example, 


• Although this assumption of exponen- 
tial dependence of the variables is inherently more 
accurate, almost identical results (not shown) 
were obtained with the approximate method by 
assuming linear histories of the variables within 

each interval ^i.e., 

The first application was to simulate the motions 
of a rocket model which was the last stage of a 
multistage rocket system. The model can be 
thought of as a cylinder about 1% feet in diameter 
and 3% feet in length having a ratio of fuel weight 
to total weight of 1/2. An angular thrust asym- 
metry of 0.001 radian in both the pitch and yaw 
planes provided a continuous disturbance to the 
motion. A separate spin motor was assumed to 
increase model spin rate during the problem from 
5 to 9 radians per second. The approximate 
solution was computed with two intervals and 
with ten intervals, each with and without jet 
damping. Problem constants and initial condi- 
tions are listed in appendix B. Results are shown 
in figures 6 to 8. 

Figure 6(a) shows the approximate 6 and \p 
histories obtained with the two-interval solution 
(no jet damping) and their comparison with the 
numerical solution. In general, this comparison 
indicates a good approximation of the numerical 
solution except for the first negative peaks of 
each curve where the approximate solution under- 
estimates the actual values. The phase difference 
between the approximate and numerical solutions 
is to be expected and is usually of little importance. 
In this respect, both ends of all intervals are 
exactly in phase, the greatest difference occurring 
halfway through each interval. The two-interval 
solution for flight-path direction is shown in figure 
6(a) along with the numerical solution. Here, 
the comparison appears not quite so good as the 
attitude solution, but satisfactory for most 
purposes. 

In order to illustrate the accuracy obtained 
with many intervals, the ten -interval solution of 
figure 6(b) was computed. Note the improve- 
ment in the X and y solutions as compared witli 
the two-interval results. 
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(a) Two-interval approximation. 

Figure G. — Approximate and numerical solutions of first sample problem. j= 0. 
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(b) Ten-interval approximation. 
Figure 6. — Concluded. 
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Thus far, problem cases have been restricted to 
zero jet damping because the program used to 
obtain the numerical solutions does not include 
jet -damping effects. However, jet -damping ef- 
fects can be quite large as illustrated in figure 7. 
This figure presents the ten-interval solution with 
and without jet damping. Note the large attenu- 
ation effect of jet damping on the low-frequency 
(precession) mode and the near absence of this 
effect on the high-frequency mode. This large 
damping effect is not unusual and can even be 
called typical. Its presence is fortunate since 
no other forces are naturally available to damp 
the motion. 

It is well documented that bodies having I x = 
I Y —J Z cannot be spin stabilized. However, the 
motions of axisymmetric bodies passing through 
this condition are not well known. For this 
reason, the second type of problem for simulation 
was selected to reveal the effects of passing through 


inertial resonance 



This problem as- 


sumed a variable I x , a constant /, and no roll 
inputs. Only the X solution was computed for 
purposes of simplicity. The problem was com- 
puted with one interval, two intervals, and three 
intervals, in all cases with no jet damping. The 
necessary constants and initial conditions are 
presented in appendix B. Results are presented 
in figure S which also shows the numerical solution 
for comparison purposes. Figure 8 illustrates the 
one-interval solution ol this problem to be inade- 
quate. The two-interval solution is much im- 
proved and reveals the trends of the numerical 
solution. However, for an accurate amplitude 
comparison, the three-interval solution is 
indicated. 

Intervals. — Solutions have been previously de- 
scribed as one interval, two interval, and so forth, 
with no explanation of why or how the number of 
intervals was selected. As previously mentioned, 
a two-interval solution (for example) means that 
the problem is computed in two intervals, usually 
so as to result in about the same percentage changes 
of the variables within each interval. Then, 
closed solutions for the first interval are obtained 
with equations (16) and (25) along with the initial 
conditions of the problem. In order to compute 
the solutions at any time of the second interval, 


however, the final values of interval one must be 
calculated and used as initial conditions for the 
second interval. 

While the number of computations increases 
with the number of intervals used, the accuracy of 
the results increases thereby. The optimum 
number of intervals, then, depends upon the 
computing facilities available, the degree of 
accuracy desired, and the total percentage change 
of variables throughout the problem. In this 
respect, the author has tentatively settled on 


using intervals in which the values of p or p~ 

do not vary more than about 15 percent for 
“accurate” results or more than about 30 percent 
for approximate results. These percentages are 
based on a limited amount of experience. Per- 
centages for the sample problems are given in 
appendix B. 

Computing time. -Computing times for sample 
problem number one of appendix B were obtained 
for both the approximate and the numerical or 
step-by-step solutions. These solutions were gen- 
erated by an IBM 7090 electronic data processing 
machine and required certain compatibility 
changes since both programs were originally set 
up for the IBM 704 electronic data processing 
machine. 

Both programs required about 2G seconds read- 
in time. Excluding read-in time, the approximate 
method (10-interval solution) had a ratio of ma- 
chine time to problem time of about 0.93 and the 
numerical method a ratio of about 4.8. Reducing 
the number of intervals used in the approximate 
solution would decrease its ratio only a small 
amount. 

Other factors involved in a computing time 
comparison are as follows: First, the numerical 
method is programed to yield output quantities 
not included in the output of the approximate 


method. It is estimated that the elimination of 
this part of the program would amount to about 
% reduction in computing time for the numerical 
method. Second, the above ratios are for defining 
the output every 0.1 second and could be reduced 
proportionally for the approximate method by 
using fewer output times. The numerical method 
would not benefit in this respect. 
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Figure 7. -Effects of jet damping in first sample problem. Ten-interval solutions. 
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(a) One-interval approximation. 

8. — Comparison of approximate and numerical solutions of second sample problem showing motions of a bod 
passing through inertial resonance = at 2 = 10. j= 0. 
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(b) Two-interval approximation. 

(c) Three-interval approximation. 
Figure 8 . — Concluded. 
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CONCLUDING REMARKS 

A method for approximating the vacuum mo- 
tions of symmetrical rigid bodies with nonconstant 
spin rates and inertias has been developed. The 
analysis includes the effects of time varying thrust 
misalinements, mass unbalance, and jet damping. 
The method was derived for bodies having equal 
moments of inertia about their pitch and yaw axes 
and is based on body pitch and yaw attitudes being 
limited to “small angle” oscillations. 

Results have been presented in the form of equa- 
tions for space-referenced Euler angles, flight-path 
angles, and earth-referenced vehicle-trajectory 
coordinates. Equations for determining maxi- 
mum wobble have been developed for certain input 
conditions. Also, equations for body-referenced 


attitude rates, angle of attack, and angle of sideslip 
are included for convenience. 

The general solutions give insight into the indi- 
vidual effects of the variables and, in many cases, 
offer a quick means for obtaining approximate 
solutions. Although the method is somewhat 
lengthy for accurate hand computation in most 
cases, it is readily programed for automatic com- 
puter solutions. 

The method has been shown to compare closely 
with numerical solutions of two sample problems. 
The sample problems also illustrated the relatively 
large effect of pitch and yaw jet damping on body 
motions. 

Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Station, Hampton, Va., April 24, 1961. 


APPENDIX A 


EXACT SOLUTION WITH NONCONSTANT SPIN RATES 


Combining equations (8), (11), and (13) results 
in the solution 

X=X 0 +e‘>»{ A'(l +«*)«'• (?) ,0g ' (1+a ° 
-A'-B'-C'-D'+W'il+at ) 2 

+C'(l+at) 3 +D'(l+aty)e i ^ l0Se U+ '" ) } 

where 

A_B.C 

X 0 e~ i<t>0 a a 2 'a' 

a+ip„Y + 

B_2C C 


A B .C 

CL CL 2 Q ft 

[ 1 +<j)( 1 -T)] <2 «+ i ' , - ) 

B 2(7 


[ 2+i («‘)( 1_ T)] (3 “ +i * >J 

c 

a? 


This result displays the spiral nature of the 
6, \p motion for nonconstant spin rates. Spin rates 
decreasing with time result in spiral motions of 
increasing magnitude and spin rates increasing 
with time act to reduce the magnitude of the 
0 , yf/ motion. 
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APPENDIX B 

SAMPLE PROBLEM INFORMATION 


PROBLEM 1 


PROBLEM 2 


Problem length, sec 20 

Mass, slugs 20— ^+f-|oo <2 

Thrust, lb... 5, 000 - 501 

Pitch and yaw inertia, slug-ft 1 2 3 25. 6 — 0. 532 

Roll inertia, slug-ft 3. 93 — 0. 09652 

Roll input moment, ft-lb 0. 572 

Pitch input moment, ft-lb_ 0. 001(5, 000 — 502) (2+0. 0252) 
Yaw input moment, ft-lb __ 0. 001(5, 000 — 502) (2 + 0. 0252) 

Thrust arm, ft 2 + 0. 0252 

Initial flight-path velocity, ft/sec 5, 000 

Initial roll angle, radians 0 

Initial pitch angle, radians 0. 05 

Initial yaw angle, radians 0. 04 

Initial roll rate, radians/sec 5 

Initial pitch rate, radians/sec 0. 015 

Initial yaw rate, radians/sec 0. 02 

Initial flight-path angle in vertical plane, radians.^ 0. 02 
Initial flight-path angle in horizontal plane, radians, 0. 01 
Two-interval solution : 

Intervals 0^2^12, 12<2^20 

Maximum change of p within the interval, 

percent 41 

lx 

Maximum change of p within the interval, 

percent 33 

Ten-interval solution: 

Intervals 2-sec intervals successively 

from 0 to 20 sec 

Maximum change of p within the interval, per- 
cent 6 

Maximum change of p — within the interval, 
percent 4. 6 


Problem length, sec 20 

Pitch and yaw inertia, slug-ft 2 25 

Roll inertia, slug-ft 2 31. 25 — 0. 69732 

Roll moment input, ft-lb 0 

Pitch moment input, ft-lb — i 

Yaw moment input, ft-lb 0 

Initial roll angle, radians 0 

Initial pitch angle, radians 0 

Initial yaw angle, radians 0 

Initial roll rate, radians/sec 5 

Initial pitch rate, radians/sec 0 

Initial yaw rate, radians/sec 0 

Initial flight-path angle, vertical plane 0 

Initial flight-path angle, horizontal plane 0 

One-interval solution : 

Maximum change of p within the interval, 

percent 0 

I x 

Maximum change of p -y within the interval, 

percent 36 

Two-interval solution: 

Intervals 0^2^12, 12<2^20 

Maximum change of p within the interval, 
percent 0 


lx 

Maximum change of p y within the interval, 


percent 20 

Three-interval solution : 

Intervals 0<2^6. 7, 6.7^2^13.3, 13.3^2^20 

Maximum change of p within the interval, 

percent 0 

Ix 

Maximum change of p -y within the interval, 


percent 
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